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Abstract 

In this work we address the question of proving the stability of elliptic 2-periodic orbits for strictly convex 
billiards. Eventhough it is part of a widely accepted belief that ellipticity implies stability, classical theorems show 
that the certainty of stability relies upon more fine conditions. We present a review of the main results and general 
theorems and describe the procedure to fullfill the supplementary conditions for strictly convex billiards. 
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• <-h 1 Introduction 



Let a be a plane, closed, regular and strictly convex curve. The billiard problem on a consists in the free motion of 
a point particle in the plane region enclosed by a, with unitary velocity and being reflected elastically at the impacts 
with the boundary. The trajectories are polygonals in the region. 

a(s,) 



4 a(s,) 




The motion is completely determined by the point of reflection at a and the direction of motion immediately after each 
reflection. For instance, the arclength parameter s, which locates the point of reflection, and the tangential component 
of the momentum p = sin 8, where 9 is the angle between the direction of motion and the normal to the boundary at 
the reflection point, describe the system. Good introductions to billiards can be found in 0], [Hj, [El, |10j - [T5| or 

M- 

The billiard model defines a map T which to each (s,p) in the annulus A = [0, L) x (—1,1), representing the pair 
impact coordinate and direction of motion, associates the next impact and direction: 

T : A — > A 

(s,p) .— » (S(s,p),P(s,p)) 
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Since the particle can travel along the same polygonal in both senses, the problem is time-reversing and the inverse 
map T -1 is well defined. 

The derivative of T at (s,p) is implicitly calculated and is given by the formulae: 

dS _ l(s,p)-R(s) cos 6(p) 
~ds ~ R(s) cos 9{P) 
dS = l(s,p) 

dp cos 9 (p) cos 9 (P) [ ' 

dP l(s,p) - R(s) cos 9(p) - R(S) cos 9{P) 

ds R(s)R(S) 
dP _ l(s,p)-R(s) cos 6{P) 
~dp ~ R{S) cos9{p) 

where S stands for S(s,p) and P for P(s,p), l(s,p) is the distance between the two consecutive impacts points ct(s) 
and a(S), R is the radius of curvature of a and cos 9(p) = y/l — p 2 is the normal component of the momentum. 

If a is a C k curve, k > 2, the billiard model gives rise to a discrete two-dimensional C k ~ 1 area preserving dynamical 
system, whose orbits are given by 

0(a,p) = {T*(8 )P ),j€Z}C A. 

A billiard has no fixed points. However, given n > 2, Birkhoff's Theorem states that T has at least two different orbits 
of period n which will be fixed points of T n . The linearization of T n at any of these fixed points, say (s,p), gives the 
linear area preserving map DT™ s p «, which has a fixed point at the origin (0, 0). According to the eigenvalues of this 

linear map, the fixed point (s,p) is classified as: hyperbolic if the eigenvalues of are p and p, G M, p ^ ±1, elliptic 
if the eigenvalues are p = e %1 and ~p, p 2 ^ 1 or parabolic if the eigenvalues are 1 or -1. 

In the hyperbolic case, the Hartman-Grobman Theorem (see, for instance, [S] or ^21) assures that, on a neighbourhood 
of the fixed point (s,p), the dynamical behaviour of T n is the same as the dynamical behaviour of DT™ s ^ on a 
neighbourhood of the origin. So, (s,p) is an unstable fixed point of T n and {(s,p), T(s,p), T n ~ 1 (s,p)} is an 
unstable periodic orbit of T. In this case, the instability of the equilibrium of the linear map DT™ s ^ implies the local 
instability of the periodic orbit for the complete map T. 

In the elliptic case, the linear map DTP p j is a rotation: the origin is surrounded by closed invariant circles and is a 
stable equilibrium. However, this beautiful behaviour may not be inherited by the map T n , as it can be seen in the 
examples on the figure bellow. For both of them, the fixed point is linearly elliptic. On the left side, the non linear 
map exhibit invariant closed curves surrounding the fixed point, which is then stable. For the non linear map on the 
right, no invariant curves can be observed and the fixed point seems to be unstable. 1 




Moreover, it is not even clear if the pictures above, obtained by numerical simulations, correspond to the true behaviour 
of the maps. In fact, Dias Carneiro and ourselves [S] have proved that any C 1 strictly convex billiard map with 

1 Even more surprising is the example given by Anosov and Katok in I of an ergodic area-preserving map of the disc \z\ < 1, with an 
elliptic fixed point at z = 0. The ergodicity implies that the fixed point is unstable. This example does not represent a billiard map and 
we don't know if there are any billiards with this property. 



2 



an elliptic 2-periodic orbit can be approached by billiard maps with a 2-periodic orbit surrounded by closed invariant 
curves, ie, with a stable orbit. We guess this result can be extended to any period. Therefore, because of natural 
numerical round off errors, one can not be sure that the simulation corresponds to the actual billiard at and not to a 
very close one. 

As a consequence, in the elliptic case, a more carefull approach is needed and higher order terms must be taken into 
account to assure the local stability of periodic orbits. A classical way to handle this problem is to use the Birkhoff 
Normal Form and Moser's Twist Theorem |14| . 

In what follows we explain how this can be performed and applied to the billiard map in the case of 2-periodic elliptic 
orbits. We have employed the software ©Maple to calculate the necessary data and all the worksheets are available 
at http:/ /www. mat. ufmg.br/~syok/papers We then apply the results to two special classes of billiards. 

Related works are Hayli, Dumont, Moulin-Ollagnier, Strelcyn p] and Moeckel The first authors studied the 

stability of periodic orbits for a family of Robnik's billiards. The last author studied the generic behaviour of the first 
Birkhoff coefficient for one-parameter families of conservative maps. 



2 Non linear analysis and the local stability of elliptic orbits 

Let T be an area preserving map with a n-periodic orbit {(0, 0), T(0, 0), T" -1 (0, 0)}. We will assume that the map 
is C k with k > 4. In the case of the billiard map, this is equivalent to assume that the curve a is at least C 5 . 

The map T n can then be expanded in Taylor Form up to order 3 in a neighbourhood of its fixed point (0, 0), 

T n (s,p) = (a 1Q s + a Q1 p + a 20 s 2 + ... + a Q3 p 3 , b 10 s + b 01 p + b 20 s 2 + ... + b Q3 p 3 ) + C(|(s,p)| 4 ). (2) 

If the fixed point is elliptic, with eigenvalues /i = cos 7 -I- isin7 and p,, by means of a complex linear area preserving 
coordinate change which diagonalizes the linear part, the map T n can be written as 

z 1 ► fi(z + C20Z 2 + ciizz + co2Z 2 + c 3a z 3 + c 2 iz 2 z + ci 2 zz 2 + cq 3 z 3 ) + 0(\z\ 4 ). (3) 

If fi J / 1, j = 1,2,3 or 4 we say that \i is non resonant and an analytic coordinate change brings the map into its 
convergent Birkhoff Normal Form 

z h-» e l ^ +T ^z + 0(\z\ 4 ) = ^z + i^ Tl z\z\ 2 + 0(\z\ 4 ). 

The first Birkoff coefficient t\ is given by 

ev ^ , sin T (?\ |2 , 2cos7-l 2 \ 

Tl = 3 (C2l) + 3 c 20 + — c 02 (4) 

cos 7 — 1 \ 2 cos 7+1 J 

where 3(c2i) stands for the imaginary part of C21. 

The calculations leading to expression (@J are standard (0, an d can also be easily performed using symbolic 

programming. 2 

By Moser's Twist Theorem if the first Birkhoff coefficient n is not zero there are T n -invariant curves surrounding the 
fixed point and therefore it is stable. We have that each point of the n-periodic orbit is contained in an open set, called 
an island, homeomorphic to a disk and invariant under T n . Each island contains T™-invariant curves surrounding the 
periodic point. So, the n-periodic orbit of T is stable. 

2 http: / /www. mat. ufmg.br/~syok/papers/3NormalForm. 
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3 Elliptic 2-periodic orbits of convex billiards 



Any closed regular strictly convex C 2 curve a has at least two diameters, characterized by points with parallel tangents 
and equal normal lines (like the axis of an ellipse) . The motion along each one of them corresponds to a 2-periodic 
trajectory for the billiard map associated to a. 

It is easy to prove that the longest of these diameters, if isolated, corresponds to a hyperbolic orbit (see, for instance 
[S| or [5]). The other(s) can be either hyperbolic, elliptic or parabolic. Let us suppose that one of them is elliptic and 
let s — and s = s% be the arclength parameters of the trajectory. As the motion occurs in the normal direction, the 
tangential component of the momentum p is zero in both of the reflection points. Then {(0,0), (si,0)} is an elliptic 
2-periodic orbit of the associated billiard map T(s,p) — (S(s,p), P(s,p)) and (0,0) is an elliptic fixed point of T 2 . 

Let i?o = -R(O) and i?i = R(si) be the radii of curvature of a at the s — and s = si and L = \\a(0) — a(si)|| be the 
length of the trajectory. Using equation Q}, the linear map DT 2 ^ = DT( Sl0 ) DT( Q ) is given by 



/ (L — Ri) (L — R ) | L(L ~ R — Ri 



V 



i?l Rq 



Rq Ri 



2(L-R )(L-R -R 1 ) 
Rq 2 Ri 



2L(L-R 1 ) 
Ri 

(L - Ri) (L - R ) , LiL-Ro-Rx) 



Ri Rq 



Ro Ri 



\ 



(5) 



and its eigenvalues are 



(L - Ri) (L - R Q ) 
Rq Ri 



1 ± 



2^L {L-Rg- g0 (L - R x ) (L - R ) 

Ro Ri 



As the trajectory is elliptic the relations L — Rq — Ri < and (L — Rq)(L — Ri) > must be fullfilled. Assuming that 
4 (L - Rq)(L — Ri) ^ R R 1 and2(i- Rq)(L - RJ ^ R R X then p? ^ 1 for j = 1,2,3,4. 

In the elliptic and non-resonant case, in order to investigate the stability of the fixed point, we can proceed and 
examine the first Birkhoff coefficient given by (@J . The complex coefficients C21, C20 and C02 in the formula depend 
on the real coefficients a^- and bij of the Taylor expansion of T 2 at the origin J2J). 

The linear coefficients Oy and fey, i + j = 1, are obviously the entries of -DT 2 , Q ^ and thus given by ©. Note that 
a io = box. As T is area preserving, a 2 — aoi^io = 1 and, as (0, 0) is elliptic, aoi^io < 0. 

These conditions were used to write down the coordinate change leading to and we got 3 : 

aio 



3(c 21 )= ^ 



C20 



l c 02| 



1 

16 
1 

16 
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ho 
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ao2 + a 2 o + hi) + 



ao2 + a 2 o - bn) + 



1 

16 
1 

16 



bw 
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aoi 



aoi ( 
610 
aoi . 
ho 



J20 



an 



an 



(6) 



which shows that n is linear on the real coefficients of third order and quadratic on the second order ones. 



In order to calculate explicitly the first Birkhoff coefficient, all is needed now are the second and third order coefficients 
of the Taylor expansion at (0,0) of T 2 {s,p) = (S(S(s,p), P(s,p)), P(S(s,p), P(s,p))) . 

3 http: / /www. mat. ufmg.br/~syok/papers/2Complex. 
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A sequence of straightforward but long computations using the Chain Rule gives those Taylor coefficients. To illustrate 
it, let us give the expression of 020 



a 20 = ^S(S(a,p),P(a,p))(0,0) 



ds 



ds 



dsdp 



2 ds ' ds 

ds 



£<M> 

OS 



— ( s l>°) 



dP 
ds 



92 S , n s 



The first derivatives of the functions S and P are given by formulae ijl^nd they depend on the function l(s,p). So, to 
calculate he second and third derivatives of S and P it is necessary to evaluate the first and second derivatives of I. 
Let l(s,S) = \\ct(S) - a(s)\\. Then l(s,p) = l(s,S(s,p)). 

By differentiating 

l 2 (s, S) = (a(S) - a{s), a(S) - a(s)) 

we have 

dl 

l(s, S) -(s, S) = - (a'(s), a(S) - a(s)) (7) 
and so, as a' is the unitary tangent vector, 

g- s (S,S)=-p. 

Analogously 

^(s,S) = P. 

These relations simply shows that —l(s, S) is the generating function of the billiard map, as a Twist map. 
Differentiating (JZJ) with respect to s and S and using that 7/ = R a" is the unitary normal vector gives 



dn. 1-p 2 ^rr~f 



ds 2 ' ' ' l(s, S) R(s) 



d 2 l ^_ ^(l- P 2 )(l-P 2 ) 



dsdS y ' ' l(s, S) 

The same reasoning gives 



d 2 l 1-P 2 VI - pa 

dS 2 ^ ' ~ l( 8 , S) ~ R(S) 

The chain rule will give, now, the first and second order derivatives of l(s,p). To evaluate them at (0, 0) and (si, 0) it 
is usefull to remember that a'(si) = —a'(0), rj(si) = —17(0) and a(si) — a(0) = Lrj(0). 

Because of the recurrent structure of the formulae, the explicit calculus of the ay and bij is suitable to be imple- 
mented as a computer program 4 . The final expression of the Taylor expansion of T 2 is also given in the worksheet 
ITaylorCoef f s. 

The second order coefficients <Hj,bij,i + j = 2 will have linear dependence on 4^(0) = R' and ^r(si) = R[ while the 

third order coefficients a,y, b^, i + j — 3 will have linear dependence on ^pr(O) = Rq and 4-r(si) = i?" and quadratic 
on the first order derivatives. So the first Birkhoff coefficient t\ will be quadratic on the first derivatives of R and 

4 http://www. mat. ufmg.br/~syok/papers/OThreeJet and ITaylorCoeffs 
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linear on the second ones. The final expression of T\ is obtained after substitution of the and bij into and then 
into Q given 

^ _ 1 Ro + Ri 1 L (L — Ri i — i?o ^ w 

8 i?o -Ri 8 L — i?o — Ri \ i — -Ro " — -Ri 1 

1 L ^ L - Rl (Ro) 2 +2^(R[f + 3R^R' 1 



8 (-L — Rq — -Ri) 2 \ L — i?o -£/ — i?i 

- H a H 1 



8 (L-R Q -R 1 ) 2 (A(L-R )(L-R 1 )-R Q R 1 ) \ L - Ro R L - R 1 R x 
The details leading to the formula are given in the worksheet 4Tau. 



4 Billiards with islands 



As remarked in section |3 a billiard on a strictly convex C 2 curve always has 2-periodic orbits and the largest one, if 
isolated, is hyperbolic. Unfortunately, one cannot assure that at least one of the others is elliptic. In fact, there are 
many examples where all 2-periodic orbits are isolated and hyperbolic (see, for instance [S] or 

On the other hand, cllipticity is an open property, in the sense that if a billiard associated to a C 2 strictly convex 
curve a has an elliptic 2-periodic orbit, then any strictly convex curve sufficiently C 2 -close to a generates a billiard 
with an elliptic 2-periodic orbit So, a large class of strictly convex billiards has elliptic 2-periodic orbits. The 
question is: are they stable? 

In what follows we present two classes of billiards (locally circular and symmetric) exhibiting stable 2-periodic orbits. 



4.1 Locally Circular Billiards 



Our first and simplest example is a 2-periodic orbit between two circles. More precisely, let a be a C 5 plane strictly 
convex closed curve parameterized by the arclength s, with the following properties: 

- there are two points located by s = and s — si such that a'(0) = —a'(si) and a(0) — a(si) = —Lff(0), where 
77(0) is the unitary normal vector at 0. 

- a is locally a circle, both near s = and s = Si, with radii Rq and R\ respectively. 

- L, R and R x verify L — R Q — Ri < 0, (L - Rq)(L - R x ) > and 4(L - R )(L - Ri) j= R Ri, 2i? i?i, which are 
open conditions on the (L, Rq, i?i)-space. 




G 



With these properties, {(0, 0), (s\, 0)} is a non-resonant elliptic 2-periodic orbit for the billiard map T associated to a. 

As a is locally circles, T is locally analytic and the first Birkhoff coefficient of the elliptic orbit can be calculated. 
Moreover, R' and R" vanish at s — and s = s\. So, 

8 (i?o Ri} ^ 

and this billiard has a stable 2-periodic orbit. 

Although extremely simple, this example shows that exchanging the curve a by the osculating circles at the impact 
points gives information about cllipticity, but not about stability, since n depends on the derivatives of the radius of 
curvature. 



4.2 Ovals with a special symmetry 

Let R be a periodic C 4 function with Fourier expansion 

R((fi) = ao + a n cos 2nip 

71 = 1 

with a n > and a > J2n=i a «' implying that R(<p) > 0, VV>. 
Let a be a curve, having R as its radius of curvature, given by 

a{<p) = {x(<p),y(<p)) - R((3) cos /3df3,£ R((3) sin f3df3 

It is a regular, closed and strictly convex C 5 curve. 

y 



OC(7t) 




oc(O) 



As R is an even function, x(— <p) = —x(ip), y(—<f ) = y(f ) and a(0)a(ir) is an axis of symmetry for a. {(0,0), (n, 0)} 
is a 2-periodic orbit for the associated billiard map. 

We have 

L = "°W -«™ -**> " 2 °°-g (2n + l)pn-l) 
J2(0) = i2(7r) = Ro = a Q + J2 a n 

n=l 
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and then 



L - R(0) - R(tt) = L - 2R = -2^< 
L - R(0) = L - R(n) = L - R a = a 



1 



-^2 an 



n=l 



(2n + l)(2n - 1) 
2 

(2n + l)(2n - 1) 



1 < 



If 



(2 + fc) a - Y, 



+ (2-fc) 



^0, fc = 0,±l,±\/2 



_(2n + l)(2n- 1) 

which are open conditions, then {(0, 0), (it, 0)} is elliptic and non-resonant. 

Let s = s(tp) be the arclength parameter for a. Choosing s(0) = and s(tt) = Si we have 



d 2 R 
ds 2 



dR 

ds 

1 d 2 R 



1 di? 



= 



Rl dip 2 



<p=0 



ip=0 

-4 ^ " 2a " < 



n=l 



ds 2 



1 rfi? 

i?0 



= 



if — 7T 

1 d 2 i? 



= -4 ?i 2 a„ < 



n=l 



and the first Birkhoff coefficient is 



So the 2-periodic orbit is stable. 



Tl 



I 

4i^ 



d 2 i? 



R Q (L-2R ) dip 2 



(0) < 0. 



In particular, this class of curves include those studied numerically by Berry in and defined by R(ip) = 1 + e cos 2cp 
with < e < 1. If e ^ |, ^ or ^-(13 — &v2), the conditions for non-resonant ellipticity are fullfilled and the 2-periodic 
orbit is stable. It is claimed in [2] that when e = | (meaning parabolicity of the 2-periodic orbit) there is neutral 
stability. There is no specific observations for the other two values of e. 

It would be interesting to investigate the behaviour of this and other examples at resonances. 
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